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SUI^BY 

A Bisthod of designing vaneless diffusers, using data given for 
simpld conical diffusers, is presented. The rate oi expans...on of 
the flow area through the diffuser is foi.ind in terms of an equiva- 
3 ent cone laid off along the logarlthaiic spiral having an angle 
equal to that for the optimum frow condition of the Impeller-. 

A family of diffusers with equivalent cone angles of 4°, 6°, ^ 
and 6®, all having the sane throat height, was designed and experi- 
mentally studied. A second series of diffusers, having throats 
equal to 62, 72, and 93 percent of the height just beyond^the 
impeller outlet, and all having the cone angle corresponding to the 
best performance in the first tests (6°), was then constructed. A 
diffuser with the 6° equivalent cone angle and a contraction ratio 
of 0.72 was found to have the highest efficiency, about 0.83 at the 
diffuser exit at tip speeds of 900 to 1200 feet per second. The 
improved performance with the forementioned design parameters was 
particularly noticeable at high loads and tip speeds. 


INTRODUCTION 

The demand for improved diffuser performance becomes increasingly 
insistent with refinement in impeller design and emphasis on high 
compressor efficiencies. As a device for achieving high efficiency 
in the centrifugal-type compressor over a wide range of air flows 
and tip speeds, the vaneless-type diffuser has certain inherent 
advantages over the vaned-type diffuser . At the impeller outlet 
the distribution of velocities, both in magnitude axid direction, 
precludes the possibility of introducing fixed vanes that will 
provide the optimum angle of attack for all the incident air at any 
one operating condition, and changes in this angle with volume flow 
and tip speed still further reduce the probability of satisfactory 
incidence. Moreover, the presence of diffuser vanes near the impeliei 
outlet may result in relative velocities equal to that of sound, with 
consequent shock losses. On the other hand, continuous diffusion oi 
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supersonic circ^uiferential-velocity ccnip oneuts apparently can be 
accomplished in the vaneless diiTuser and cnanges in the flow angle, 
in the absence ol* vanes, appeirantly will produce no serious increases 
in diffuser losses. 

The design of vaneless diffusers involves two closely related 
problems. In order that the diffusion may be efficiently accomplished, 
the nonuniform flow at the impeller outlet should be converted into 
a steady flow with a uniform velocity profile before a serious 
attempt is made to reduce the axial and radial components of velocity. 
This transformation is a function of the transition section, the 
passage from the impeller outlet to the narrowest part (throat) of 
the diffuser passage. The function of the diffuser proper is to 
convert the dynamic energy of the advancing flow into pressure as 
rapidly as possible without incurring flow separation at the most 
adverse operating condition. Because diffusion of the tangential 
component' is essentially a function of radial distance only, any 
acceleration of the diffusion process must be effected by increasing 
the rate of diffusion of the radial and axial components of velocity, 
that is, by increasing the rate of expansion of the flow area. 

An efficient diffuser should perform both of these operations 
with the smallest possible losses and at the same time, because of 
the penalty on size, in the sDiallest possible frontal area. Inas- 
much as detailed information concerning flow conditions in the transi- 
tion section is lacking, only rather general design criteria can be 
set up; for example, gradual changes in direction of the walls and 
gradual changes in flow area, particularly from small to large areas. 

In the absence of satisfactory theoretical criteria for the determina- 
tion of the optimum rate of diffusion in a vaneless diffuser, the 
optimum rate of area expansion, as a first approximation, can be 
made the same for the vaneless diffuser as the optimum rate for a 
simple conical diffuser. Experiments with such diffusers (refer- 
ence l) have showTL that maximum efficiency is associated with 
divergence angles of 5^ to 8^. 

The analytical problem is thei'’efore; (a) to determine a transi- 
tion profile that, conforming to the preceding criteria, blends 
smoothly with both impeller outlet and diffuser entrance, and (b) to 
fix such values for the distance between the two plates of the 
diffuser proper that the rate of area expansion along a prescribed 
flow path shall be the same as the rate of area expansion of an 
equivalent cone having the optimum divergence angle, vrnon equations 
incorporating the conditions on transition and diffuser sections are 
set up, a variable parameter may remain undetermined in each case. 

The optimum value of the appropriate parameter for each section may 
then be experimentally determined. 
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The iTiethod used In the design of the diffusers and the results, 
in terms of compressor* performance, of vai’ying the parameters asso- 
ciated with the transition section and the diffuser proper are 
presented. Five 34-inch-diam.efcer vaneless diffusers were studied. 
Three diffusers were designed with the same transition-section 
parameter hut equivalent cone angles of 4°, 6°, and 8 . The transi- 
tion section of the diffuser showing the best performance was then 
so modified as to give first a wider, then a narrower, throat than 
the original, and the effect of these modifications was experimen- 
tally investigated. 

SYMBOLS 

The following symbols are used: 


c 

constant of integration 

"p 

specific heat at constant pressure, foot-pounds per 
pound °F 

D 

hydraulic diameter, inches 

f 

friction coefficient (ratio of shear stress to dynamic 
pressure) 

g 

acceleration of gravity, 32.14 feet per second per 
second 

h 

passage height perpendicular to rear shroud, inches 

K 

constant dependent on ag 

1 

distance measured along flow path, inches 

M 

angular momentum per unit mass (rU^), feet per second x 
inches 

P 

total pressure, inches of mercury absolute 

P2/fl 

compressor pressure ratio, ratio of outlet total pres- 
sure to inlet total pressure 

P 

static pressure, inches of mercury absolute 


corrected volume flow, cubic feet per minute 
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,/Ve V specific capacity; cubic feet per mirwte per square foot 
’ (whoi'e Dp is iiapelier-outlet diaiueter) 




R 

gas constant, 55.4 

r 

radial distance from axis, inches 

^c 

radius of equivalent cone section, inches 

i’t 

radius of curvature for transition-section front shroud, 
inches 

T 

total temperature, 

t 

static temperature, °R 

U 

absolute air velocity, feet per second 

Ua 

component of air velocity in axial direction, feet per 
second (fig. l) 

Um 

meridional component of air velocity (normal to U,(; 
and parallel to rear sliroud of diffuser), feet per 
second (fig. l) 

Ur 

component of air velocity in radial direction, feet 
per second (fig. l) 


U 0 component of air velocity tangential to circle of 

radius r, feet per second (fig. l) 



✓ ... - W ' 

V 

impeller tip speed, feet per second 

w 

mass air flow, pounds per second 

a 

flow angle (angle between U and Ug), degrees (fig. 1) 

3 

angle betwee.n diffuser rear shi’oud and radial direction, 
degrees (fig. l) 

7 

ratio of specific heats 

6 

half -angle of equivalent cone along design path, degrees 

€ 

angle between front and rear shrouds at impeller outlet, 
degress 
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lad adiabatic efficiency 

9 ratio of actual inlet stagnation tenrperature to MCA 

standard sea-level tejiperature 

^ contraction ratio for transition section (hg/liQ) 

P air density, pounds per cubic foot 

ang] e teiveen froxit and rear stirouds at diffuser 
entrance (after ti'ansition section), degrees 

Subscripts : 

9 entrance to diffuser section (transition-section exit) 

o impeller outlet (transition-section entrance) 


DESIGN THEOEY 

The two phases of diffuser operation, transformation of the 
chaotic flow at the impeller outlet into a uniform stream and 
conversion of the kinetic <?nergy of the radial and axial components 
of velocity into pressure, may conveniently be separately considered, 
although neither component of the diffuser can be designed independ- 
ently of the other. The diffuser proper will be considered first; 
then the requirements for a smooth transition section from impeller 
to diffuser will be established. 

Di ffuser section. - V/ith respect to the diffuser proper, the 
problem as defined is the determination of a diffuser profile that 
will provide a rate of increase of area along a predetermined flow 
path the same as that of a conical diffuser of divergence angle 26. 
Because of compressibility, designing for only one operating condi- 
tion is possible. At this design condition, the air leaves the 
transition section at some angle that is assumed to remain 

constant throughout the diffuser. The design path is therefore a 
logarithmic spiral, or more precisely, because of an axial component 
of flow in this instance, a helix on a cyj-inder whose directrix 
curve is a logarithmic spiral (fig. 2), 'A straight line equal in 
length to the arc length of the helix forms the axis of the equiva- 
lent cone. 

At any distance I along the design path 



2nrh sin ag 


( 1 ) 
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where 2Ttrh sin is the area of the diffuser passage normal to 
the flow direction and Jtrc"^ is the area of the coiresponding sec- 
tion of the eq.ui’''alent cone. The quantities r, h, and r^ are 
all functions of I . 

When equation (l) is differentiated with respect to Z, 


d(rh) *^^0 

sin ~ = r^ 


( 2 ) 


But 


dr^/dZ = tan 6 


and 


d(rh) 

dZ 


d(rh) dr d(rh) 


dr dZ 


clr 


sin ttp cos P 


If these substitutions are made in equation (2), the value of 
introduced from equation (l), and t'ne equation Is rearranged, 


dr 


(rh) 


^\/2 tan 6 


lin^/^ cLq cos 3 


= K 


sin" 

where K is adjustable through its dependence on ag. 


rc IS 


(3) 


If equation (3) is rewritten in the form 


and integrated, 


Kdr 


d(rh) 

^i/rh 


Kr = 2''|rh - C (3a) 

where C, a constant of integration, depends upon the profile and 
the height of the transition section, which fix hg, rg, and ag. 
It is determined by the conditions that the diffuser entrance shall 
have the same height as the transition-section exit, and that the 
two sections, at their Junction, shall have common tangents at both 
front and rear shrouds . 
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When equation (3a) Is 


h 


solved f or h 



(4) 


Equation (A) gives, in terms of constants that m.ust be adjusted to 
the transition section, the value that h must have at all radii in 
order that the rate of area expansion along the design path shall be 
the same as that of a ccae of half -angle 5 . 


Transition section. - Pi’imarily because of the extreme com- 
plexity of the flow conditions at the impeller outlet, it is diffi- 
cult to impose on the transition section any more definite require- 
ment than that it shall have gradual changes in both the direction 
of the vails and the flov area. For simplicity this study is limited 
to a section whose front vail is a circular arc, whereas the rear 
wall, after turning gradually tiirough an angle of 5° in a radial 
direction, continues as a straight line. The 5° bend was neces- 
sitated by the dimensions of the experimental setup. The require- 
ment for smooth transition from impeller to diffuser Implies five 
geometrical conditions: 


(l) The relation between passage heights at entrance and exit 
of the transition section may bo expressed 

h0 = AIiq (S) 

where A is hero regarded as an arbitrary constant whose optimum 
value is to be determined experimentally. 


(2) The fundamental equation for the determination of h 
(equation (4)) becomes, at the exit from the transition section 


h 


e 


^^^e ICC C^ 
4 ^ 2 ^ 4r0 


(4a) 


(3, 4) Fi’om the geometry of the configuration shown in fig- 
ure 3, 


re - ro = r^(sin £ + sin q?) cos 3 


(6) 


and 


h 


r^ cos TO = Iiq + r^^ cos € 


(7) 
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(5) The req.ulreinent that the diffuser and transition sections 
have the saire slope at the diffuser entrance yields 

ten^_dh C_2_ /g\ 

cos p “ dr ~ 4 " 4^ 2 

The four physical conditions to be met are tho continuity equa- 
tion 


2jtrhUj^ 


(9) 


the adiabatic relation (appendix A) ^ 


P 




x2.53 





the equation for the loss of total piressure (appendix B) , 


( 10 ) 


dP 


^ 2 
f CPC a sec 

g ■ liI?T ■ 


a dr_P 


( 11 ) 


and the angular-momentum equation (reference 2), which in the present 
notation becomes 


^ _ f dr 

M he cos p 


( 12 ) 


where the magnitude of f is one-fourth the magnitivie of the friction 
factor of reference 2. The factor cos p appears here in the 
denominator because the element of design path of reference 2 is in 
a plane normal to the axis. 


DESIGN PEOCEDUEE 


The geometry of the impe?.ler outlet ^ as well as the mean fiow 
conditions and the state of the air at that location, is assumed to 
be known; that is, , the quantities ig, hg, W, T^ = 13 ^ Pqj ^e,o^ 
and V are given, ' From these quantities, the meridional component 
of the air velocity ^ and the ai^gle a^, where 


U. 

tan ag = 




m,o 
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can be founcl as also can the density p^. The successive steps in 
the assign are then as fo±lorfs; 

1. The value of hg vas detenained froai the relation hg = Mig. 
The initial value of A vas necessarily arbitrary. For the condi- 
tions of this investigation, the cptimum value of X subsequently 
was experimentally found to be approximately 0.7. 

2. Inasmuch as a direct solution for the uiaknown variables 
leads to excessively complicated eq.uations, a trial-and-error 
approach is more practicable. The angle cp is necessarily small 
because of the small values of 6. As a first approximation to the 
solution, cp may thei’efore be assumed equal to 0. Equations (6) 
and (?) then take the form 


where 

G 

From the preceding equations 

r„ - r^. = 7 —^ sin e cos 3 

® ° 1 - cos € 

This value of r^ - was retained and is later used to find the 
coi-rected value of cp. 

3. Equation (12) was integrated from entrance to exit of the 
transition section to determine the angular momentum of the flow at 
the transition exit. This integration provides the relation 

log Mg = - f -—-^(rg - rg) -f log (13) 

From continuity consideretions and the geonetry of the flow 


= G sin € cos 3 


^o ” ^e 
1 - cos € 


CSC 

h" 


23T 


The magnitude of rUg, and of rU, decreases slowly through the 
transition section whereas that of p increases slightly. The 
magnitude of (esc a) /h was therefore assumed to be constant through 
the section and the initial value was used. 
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4. The circumferentlaD - /eloclty component at the exit of the 
transition section vas fount from ohe equation 

^e 

^e.e - 

' ^ e 


5. When the va'ue of M from equation (13), 


. f - r^) 

M = Mq e ^ (14) 

is substituted in equation (ll) and the equation integrated between 
the limits Pq and Pg, the value for Pg is given in the form 


log Pg = 


f CSC a sec^ a 

g h KT cos 3 


J 


_ CSC a/ V 

e - l-o) 



2 \ 


_u_ 

2gCpT 




dr + log Pg 


(15) 


In equation (15) and the following equation, an approximate value 
of U is needed, A sufficient!;^ accurate value can be obtained 
from the equation + Ug'" by using the value of Ujj> at the 

transition-section entrance. 


6. The density Pg at the exit of the transition section was 
found from the relation 


.322 


Pe = 



\ 

2gCpT 




7. Values of ,cp and a can then be calculated 
following equations, which are developed in appendix 


from the 
C: 


tan cp = B 4 G cot ag 


(16) 


B = 2 tan 


1 

2 


c 


2 hg cos (3 



14 4 W co s B 

■ :tPeMg(rg -TJ 


where 
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and 


144 W cot ttg 


tan^ 8 tan n 




+ 



tan 6 


2 tan^ 8 

3in^ ttg cos^ 0 


4 tan cp 

COB 3 


(17) 


A trial-and-error inothod was employed In the solution of these equa- 
tions. A value for cp was assumed and equation (16) was solved 
for . 'J'hese values were then inserted in equation (17) and, by 
inspect] on, a second approximation to cp was chosen, such as to 
reduce the existing inequality, and the procedure was repeated. 

Tliree or four trials were usually sufficient to give an identity 
for equation (17). 

8. The constants K and C, which determine the profile of 
the passage, were then obtained from equations (3) and (8), res- 
pectively: 


9. As a check, the channel height was then calculated from the 
two equations 


tan 8 



and 



The value of r^ was found from equation (6) 


(sin £ + sin c?) cos 3 


hg = + r.j.'^cos c - cos q?) 


■"e 4 2 4rg 



(4a) 
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DS3TGWS TESTED 

The design theory laeves andeterrjined the optimuia v^.'.lues of 
both the contraction rr.tio A and the e'iui'’’aleiit cone angle 6 . 

It was therefore desirable that the e.xperirDenta.1 program investigate 
the effect of vaiying each of these parameters. Three vaneiess 
diffusers of o4-inch diameter were bui].t and tested. For each of 
these diffusej’s, the value of A was 0.72 and that of a was 26°, 
corresponding to the best impeller performance. The equivalent cone 
angles of the three diffusers were 4°, 6°, and 6° (fig. 4). The 
diffuser that gave the best performance in these tests (6° equiva- 
lent cone) was then so modified as to have values of A, first, of 
0.95 and then of 0.62 (fig= 5). 

The same rear diffuser wall was used in all the diffusers. The 
magnitude of the angle 3 (25°) was dictated by the dimensions of 

the test rig, the corrosponding outlet angle of the Impeller being 
about 30°. Even then, this 25° angie was slightly large and the 
rear shreud had to be curved foriiard for the last 4 inches. The 
effect of this curvature near the diffuser exit appeared to be 
negligible. 

In the diffuser for which A = 0.93, the transition section was 
very short; the 6° diffuser section started about 0.7 inch, radially, 
beyond the impeller outlet. (See- fig. 5.) The transition section 
having A =0.62 necessitated abandoning the design condition that 
the impeller front shroud and the wall of the transition section have 
a common tangent at the transition entrance. The entrance passage was 
so proportioned that the resultant air velocity was constant for about 
the first 3/4 inch in the transition section, which resulted in a 
relatively abrupt change iii wall direction (fig. 5), 

Tables I to III show details of the profiles of the different 
diffusers . 


APPARATUS AWD TESTS 

A variable -component test rig (reference 5) was employed in 
this investigation. The mixed-flow impeller used is shown in fig- 
ure 6 and a shetch of the impeller and the diffuser, in figure 7. 
Standard instrumentation was used where applicable. In addition, 
total-pressure tubes, which provided data ior direct calculation of 
the efficiency of the impeller-diffuser combination, were installed 
at the diffuser exit (except in the 4° diffuser) . A single tangential 
collector outlet was used instead of two radial collector outlets, but 
previous investigations showed no appreciable difference in the 
instrument readings for the two configurations. 
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Each diffuser was investigated at actual tip 'speeds of 900,^ 
1100, 1200, and 1300 feet per second (except the 4 diffuser, which 
was not run at 1300 ft/sec), and the data were subsequently corrected 
to standard atmospheric conditions. For each tip. speed, the flow 
range was from open throttle to surge. 


RESULTS AND DISCUSSION 

In the comparison of the performance of diffusers with dif- 
ferent values of contraction ratio X and equivalent cone angle 5, 
efficiency at the diffuser exit rather than over-all efficiency was 
chosen as the most direct, practicable measure of the diffuser 
ciency. Although this procedure involves the tacit assumption that 
no change occurs in impeller performance with the changes in diffuser 
variables, measurements of total pressure at the impeller outlet, 
which would allow the isolation of diffuser performance, are, in 
general, rather unreliable. 

Efficiencies at the diffuser exit for the three diffusers 
having X = 0.72 and equivalent cone angles of 4 , 6 , and 8 , 
are shown in figure 8. The adiabatic efficiency was found by the 
procedure recommended for determining over-all efficiency (refer- 
ence 4) except that the arithmetic mean of the total pressures at 
four stations across the diffuser exit was substituted for the 
corresponding quantity at the outlet measuring station, lor the 
present purpose, efficiency is plotted agaiiist load coefficient, 
inasmuch as this parameter, being related to the flow angle a , 
facilitates comparison of efficiencies at approximately equivalent 
operating conditions. 

Inasmuch as total-pressure tubes were not installed at the 
diffuser exit of the 4° diffuser, the impeller-diffuser efficiencies 
were computed from; a collector correlation that gave the static 
pressure at the diffuser exit and the calculated exit velocities. 

The curves of figure 8 indicate relatively high performance 
of the 6° diffuser at high loads. The peak efficiency of the 6 
diffuser at high tip speeds is 4 to 5 points above that for either 
of the other two diffusers. In general, the 4® diffuser has a 
definite advantage at low loads. For the 4° diffuser, a decrease 
in load coefficient (and so of %) below the value for maximum 
efficiency seems to be associated with smaller losses in efficiency 
than for the 6° or 8° diffuser. At values of load coefficient less 
than the design value, the departure from optimum flow conditions 
in the 4° is leas than in the 6° diffuser, as the tendency to 
increase a towards its design value is greater for the smaller 
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equivalent cone angle. The peak adiabatic efficiencies at the 
diffuser exit for the three equivalent cone angles and the medium 
throat are given in the following table: 


Peak adiabatic efficiency 


Diffuser 


Actual tip speed 
(ft/sec) 



900 

1100 

1200 

1300 

40 

0.83 

0.82 

0.78 


6° 

.84 

.64 

.83 

0.79 

8° 

.85 i 

.80 

.77 

.75 


For the impeller investigated, an equivalent cone angle of approxi- 
mately 6° gives the highest peak efficiency at high tip speeds and 
the highest efficiency at the high values of load coefficient. 

Diffuser-exit efficiencies for the 6° diffuser with the three 
values for X (0.62, 0.72, 0.93) are shown in figure 9. The medium 
throat (X = 0.72) shows appreciabi;^'- higher efficiencies at medium 
and high loads at all tip speeds than either of the others. In this 
connection it may be recalled that the narrow throat made it impos- 
sible to retain at the front shroud the common tangent to the 
Impeller outlet and transition entrance, and deterioration of the 
flow at high load coefficient might be anticipated. At low load 
coefficient, the narrow throat quite consistently shows higher effi- 
ciencies than the medium or wide throats at the same value of load 
coefficient. 

The following table shows the peak efficiencies at the diffuser 
exit for the three values of X used in this investigation with the 
diffuser having a 6° equivalent cone expansion angle: 


Peak adiabatic efficiency 


Actual tip speed 
. (ft/sec) 



900 

1100 

1200 

1300 

0.62 

0.81 

0.78 

0.77 

0.71 

.72 

.84 

.84 

.83 

.79 

.93 

.81 

.81 

.78 

.75 
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Judged on the hasis of peak efficiency and efficiency at medium and 
high load coefficient, the transition section having a value of X 
of 0.72 gives the test results of the three in the present setup. 

The over -all performance of the compressor using the 6° dif- 
fuser with X = 0.72 is shown in figure 10. Pressure ratios are 
plotted against corrected volume flow in accordance with recom- 
mended procedure (reference 5) with adiabatic efficiencies shown 
as contour lines. Figure 11 compares the over-all performance of 
the compressor when the 6°, medium-throat (X = 0.72), vane lees 
diffuser was used with the performance of the same impeller with 
the manufacturer's vaned diffuser, at actual tip speeds of 800 and 
1100 feet per second. Pressures at the standard measuring station 
in the outlet pipe were used in computing these efficiencies. Data 
for the vaned diffuser were obtained on the same rig and by the 
same techniques as those for the investigation of the vaneless 
diffusers. The wider operating range of the vaneless diffuser at 
reasonable efficiencies is apparent at both tip speeds, and peak 
efficiencies are higher with the vaneless diffuser by about 
6 points and occur at somewhat higher load coefficients. 

A more extensive investigation using different types of 
impeller would be required to determine whether these results are 
generally valid. In particular, the observed variations in per- 
formance with the true flow angle a Indicate that the optimum 
cone angle in the transition section would probably be different 
for an impeller having a different flow angle at the design point. 


SUMMARY OF RESULTS 

A diffuser-design method was investigated with five vaneless 
diffusers in combination with a mixed-flow impeller in a variable- 
component test rig. The five impellers Incorporated three dif- 
ferent throat heights and three different rates of area expansion 
(equivalent cone angles) . The design method, which left throat 
height and rate of area expansion as parameters to be experimentally 
determined, resulted in compressor efficiencies equal to er better 
than those shown by the vaned diffuser designed for this impeller 
and retained the characteristically wide operating range of the vane- 
less diffuser. The diffuser with the 6° equivalent cone angle and a 
throat-inlet ratio of 0.72 was found to have the highest efficiency, 
about 0.83 at the diffuser exit at actual tip speeds of 900 to 
1200 feet per second. Throughout mosi of the useful operating range 
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the 6° equivalent cone angle gave peak efficiencies from 4 to 5 points 
higher than either the 4° or 8° equivalent cone. angle. The medium 
throat gave higher efficiencies hy 3 to 8 points than either of the 
others at the high tip speeds and load coefficients. 


Flight Propulsion Research Laboratory, 

National Advisory Committee for Aeronautics, 
Cleveland, Ohio, June 27, 1947. 
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APPENDIX A 


DERIVATION OF EQUATION FOR AIR DENSITY 
From the gas law, • ■ . 

0 - i (18) 

the relation between total and static pressure, 



(19) 


and the relation between total and static temperature,' 



an expression for density can be obtained in terms of total pressure, 
temperature, and velocity. Substitution in equation (IS) of the 
value of p from equation (19) arid that of t from equation (20) 
yields 




^ V2.53 


2gCpT 
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APPENDIX B 


DEEIYATION OF EQUATION FOE TOTAL-PEESSUEE LOSS 

The pressure loss required to overcome frictional resistance 
can be found from an extended form of Bernoulli's equation containing 
a friction term in addition to the customary terms for pressure and 
kinetic energy. 

The equation, developed in reference 6 in slightly different 
notation, is 


d_£ ^ 4fU^dZ 

P g 2gD 

From the general gas. law p = p/Et, and substitution of this value 
for p yields 

Et ^ + ~ = - 

P g 2gD 

The relation between total and static pressures 





J - -i- 


2gCpT^ 


when differentiated logarithmically yields 


^ M = 

P P 




If this equation is multiplied tiirough by Et, if the value c /e 

is substituted for 7/(7”l); and if the relation between total^and 
static temperatures, 


t 

T 


1 



is introduced, the equation becomes 
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RT ^ M = pt ^ 

P g P 

Bernoulli's equation can thus be written in terias of total pressure 


Btf = 


4fU^d I 
2sD, 


For the vaneless diffuser 


dZ = dr CSC a/cos 3 


^ sec^ a 
' ......Jr. 


t = T . 1 - 




2gCpT^ 


and 


D = 2h 

Substitution of these values gives 


dP f CSC g sec^ a 

P " ' g hKT 


M^dr 
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APPENDIX C 


DERIVATION OF EQUATIONS FOR cp AND a 

© 

The two equations in cp and can be found as follows : 


rg - rQ = r.(.(sin e + sin cp) 

cos 3 

(6) 

hg - hg = ri,(cos cp - cos 

€) 

(7a) 

Elimination of r... results in 




^e ■ 'o (sin e + sin cp) cos 3 
h'e - ho cos cp.- cos e 

Solving for Nq - h© yields 


“ ^e 


^e ~ ^0 cos cp - cos e 
cos 3 sin cp + sin c 


cos 3 


- 2 sin i(cp + e ) sin i(,q5 - c) 


2 sin c? + € ) cos P - ^ ) 


^e - ^o 
cos 3 


tan 


£ 2 _ 


_ j, tan !• € - tan ^ cp 


■ o 


cos 3 11 

1 + tan j e tan ^ cp 


For convenience, let 


A = 1 


+ tan — e 



tan 




(ho - he) A cos 3 


r • r 


then 
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Because cp is a small angle, 


and 


tan — cp 


1 tan cp 

2 


tan cp = 2 tan 


2i'iQ A cos 3 2hg A cos 3 



From equation (9) 


... 144W 144W 


and 


^eUe,e ^ Me 

Substitution of these quantities in the equation for tan cp results 
in 


tan cp = 2 tan g ^ " 


2hj5 A cos 3 

r “ r 
e o 


144W A cos 3 

- 


JtMgPe 


Ir 


cot a_ 


e 


In all actual cases, A will be very nearly equal to unity. For 

example, if € = 20° and cp = i° then ^ tan tan c? is 0.0015 so 

that A differs from 1 by only 0.15 percent. After using A = 1 as 
a first approximation, a more precise value can be foun-i, if desired, 
after cp is known but this refinement is in general unnecessary. 


If A 


is taken equal to 1, the equation for 
tan cp = B + G cot ag 


tan cp can be written 
. • (16) 


B = 
G = 


2 tan 



2hg cos 3 



144 W cos 3 
e* e ^e ~ ^o^ 


where 
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The quantities B and G can he found and give one of the two eqaa 
tions in cp and ■ 

The second equation Involving cp and is obtained from 

equation (4a) . Substituting in it the values of h from equa- 
tion (9), of K from equation (3), and of C from equation (S) 
yields 


144W 

2jtMp tan <Xq. 


2r„ tan^ 6 

O 

4 sin^ dg cos^ p 


+ 


2 sin' 


V 2 tan 5 


cos 3 


2 tan^ 6 

3 2 
sin ttg cos 3 


4 tan p 
cos 3 


2 tan^ 5 


^^e \sin^ ttg co.s^ 3 


4 tan cp 
cos 3 ; 


and, after division by rQ and rearrangement, 


144W cot a, 


e 


tan^ 6 


2rtMp r. 


tan cp 

sin^ ttQ "coa^ 3 ^ 


tan 5 


2 tan*^ 6 


i 


2 sin'^ ttg cos 3 


sin^ a. 


cos^ 3 


4 tan cp , 
cos 3 


(17) 
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TABLE I - PASSAGE HEIGHTS THROUGH DIFFUSERS WITH 


DIFFEREM' EQUIVALSHT COKE AI'fGLSS 


Radius 
to rear 
shroud 
(in.) 

Passage height, h 
(in.) 

Equivalent cone 
(deg) 

angle, 5 

4 

6 

8 

5.350 

0.606 

0.806 

0.806 

5.519 

.782 

.782 

i .782 

6.000 

.699 

.699 j 

I .699 

6.290 

.632 

.632 

.632 

6.700 

.574 

.579 

.584 

7.000 

.564 

.576 

..588 

7.500 

.549 



8.000 

.536 

.571 

.605 

9.000 

.517 

.570 

.624 

10.000 

.503 

.572 

.646 

11.000 

.492 

.577 

.666 

12.000 

.485 

.584 

.692 

13.000 

.480 

.593 

.717 

13 . 500 

.478 

.598 

j .729 

17.000 

.474 

.636 ! 

1 .820 


National Advisory Coamittee 
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TABLE II - RADII THROUGH DIFFUSERS WITH DIFFERENl’ 
EQUIVAI.ENT CONE ANGLES 


Radius 
to rear 

Radius 
to front 
shroud 
(in.) 

Mean 
radius 
(in. ) 

Radius 
to front 
shroud 
(in.) 

1 

Mean I 
radius 
(in.) ! 

1 Radius 
[ to front 
i shroud 
1 (in.) 

Mean 

radius 

(in.) 

shroud 







(in.) 


Equivalent cone angle, 6 





(defi) , 




4 

6 

8 

5 . 350- 

5.762 

5.556 

5.762 

5 . .556 

5. ’762 

5.556 

5.519 

5.907 , 

5.713 

5.907 

5.713 

5 . 907 

5.713 

6.000 

6.313 

6.156 

6.313 

6 . 156 

6.313 

6.L56 

6.290 

6.557. 

6.423 

6.557 

6.423 

6 . 557 

6.423 

6.700 

6.943 

6.822 

6.945 

6.822 

6.947 

6.823 

7.000 

7.238 

7 . 119 

7 . 244 

7.122 

7.249 

7.125 

7.500 

7 . 732 

7 R"' 6 





8.000 

8.227. 

8 . 113 

8.241 

8.120 

8.256 

8.128 

9.000 

9.218 

9.109 

9.241 

9.120 

9.264 

9 . 132 

10.000 

10.212 

10.106 

10.242 

10.121 

10.272 

10.136 

11.000 

11.208 

11.104 

11.244 

11 . 122 

11.282 

11.141 

12.000 

12.205 

12 . 102 

12.247 

12.123 

12.292 

12.146 

13.000 

13.203 

13.101 

13.251 

13.125 

13.303 

13.151 

13.500 

13.702 

13 .601 

13 . 753 

13.626 

13.808 

13.654 
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TABLE III - PASSAGE HEIGHTS AND RADII THROUGH DIFFUSERS 
WITH DIFFERENT THROAT- INLET PvATIOS 


Radius 
to rear 

Radius 
to front 
shroud 
(in.) ! 

Passage 

height 

h 

. (in.) 

Radius 1 

to front 
shroud 
(in.) 

Passage 

height 

h 

Radius 
to front 
siii'oud 
(in.) 

Passage 

height 

h 

(in. ) 

shroud 







(in. ) 


Contraction ratio, a 




0.62 

0.72 

0.93 

5.519 

5.907 

0.782 

5.907 

0.782 

5.762 

0.782 

6.000 

6.223 

.522 

6.313 

.699 

6.322 

.763 

6.290 

6.507 

.500 

6.557 

.632 

6.608 . 

.753 

6.700 

6.911 

.498 

6.945 

.579 

7.013 ■ 

,747 

7.000 

■7.210 

.497 

7.244 

. 576 

7 , 310 

. 735 

7.500 

7.709 

.496 



7.806 

.725 

8.000 

8.210 

,496 

8.241 

.571 

8.303 

.717 

9.000 

9.211 

.500 

9.211 

.570 

9.299 

.708 

10.000 

10.213 

.506 

10.242 

.572 

10.297 

.703 . 

11.000 

11.217 

.513 

11.244 

.577 

11.297 

,702 

12.000 

12.221 

.523 

12.247 

. .584 

12.298 

.704 

13.000 

13.225 

.533 

13.251 

.593 

13.300 

.709 

13.500 

13.728 

.539 

13 . 753 

.598 

13 . 800 

.712 

1 7 Hi 
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Fig. 



Figure I. - Velocity relations in vaneless mixed-flow diffuser. 
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Figure 2. - Vaneless mixed-flow diffuser showing constant angle of 
I ogar i t hm i-c-sp i ra I flow path and height of diffuser passage. 
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Fig. 3 



Figure 3. - Layout of diffuser entrance and transition section. 
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Figure 4, - Section through diffuser passage showing variations in di' 
vergence angle. (Dimensions in in. See table I.) 
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Figure 5. - Section through diffuser passage showing variations in 
throat size, (Dimensions in in. See tables H and ill.) 
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Fig. 6 
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( b) Bottom vi ew. 

Figure 6. - Mixed-flow impeller used In diffuser investigation. 






Figure 7 


Assemblj^ of mixed-flow impeller and vaneless diffuser 
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Adiabatic efficiency at diffuser exit 
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Figure 8. - Adiabatic efficiency at diffuser exit of vaneless mixed-flow 
diffuser with contraction ratio of 0.72 and three different equivalent 
cone ang I es. 
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Figure 9. - Adiabatic efficiency at diffuser exit of vaneless mixed-flow 
diffuser with equivalent cone angle of 6° ana three different values 
of contraction ratio. 


Pressure ratio, P^/P 


Fig. 10 
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Figure )0^ - Pe r fo rmance of compressor with vaneless mixed-flow diffuser 
with equivalent cone angle of 6^ and contraction ratio of 0.72. 
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8= 6°, X = 0.72 

Impeller-manufacturer's 
vaned diffuser 

Curve from figure 9 
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Load coefficient, cu ft/revolution 


II. - Performance curves for compressor with vaneless mixed-flow diffuser and with impeller- 

manufacturer's vaned diffuser. 
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